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A DIFFUSION  APPROXIMATION  MODEL 
FOR  A COMMUNICATION  SYSTEM 
ALLOWING  MESSAGE  INTERFERENCE 

Donald  P.  Gaver 
Naval  Postgraduate  School 
John  P.  Lehoczky 
Carnegie-Mellon  University 

1 . Introduction  and  Problem  Statement 

We  study  the  operating  characteristics  of  an  element  of  a complex  commun- 
ication system,  the  element  consisting  of  c channels  which  service  an 
arrival  stream  of  messages.  When  a message  arrives,  it  selects  a channel  at 
random  and  initiates  a transmission  (service)  time  of  random  duration.  If  by 
chance  the  channel  selected  is  already  occupied,  i.e.  is  in  the  process  of 
transmission,  both  messages  may  be  "destroyed,"  i.e.  terminated  before  com- 
pletion, and  the  channel  reverts  to  an  empty  or  open  condition.  The  trans- 
mitters of  the  messages  are  capable  of  detecting  the  event  of  destruction, 
and  following  such  an  event  go  into  a re-try  or  re-transmission  population, 
from  which  they  later  make  attempts  to  occupy  a channel  and  eventually  com- 
plete message  transmission.  We  develop  an  approximate  probability  model  to 
describe  the  performance  of  such  a system.  The  approximation  tends  to  become 
exact  as  the  number  of  channels,  c , becomes  large  (c-*»)  , but  numerical 
studies  indicate  that  it  may  be  quite  adequate  for  c near  ten. 

The  above  problem  bears  close  resemblance  to  a problem  of  packet  switching 
of  data  on  the  ARPANET  and,  very  likely,  on  other  satellite  communications 
systems,  as  discussed  in  Kleinrock  (1976),  p.  362  ff.  We  choose  to  represent 
the  various  interacting  populations  of  messages,  i.e.  those  in  process,  and 
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those  in  the  re-try  population,  by  means  of  stochastic  differential  equa- 
tions as  was  done  by  Gaver  and  Lehoczky  (1976)  for  a related  problem. 


2.  Mathematical  Formulations:  Model  1 


Assume  that  messages  arrive  at  a system  of  c channels  according  to  a 
non-homogeneous  Poisson  process  of  rate  cA(t)  ; that  is,  the  probability 
that  a message  arrives  in  (t,  t+dt)  is  cA(t)dt  + o(dt)  . Each  message 
then  selects  a channel  at  random;  any  particular  channel,  whether  occupied  or 
not,  is  selected  with  probability  c"^  . If  the  channel  is  unoccupied  the 
message  begins  transmission;  the  probability  that  it  completes  in  time 
(t,  t+dt)  is  p(t)dt  + o(dt)  . Note  that  if  y(t)  = p , constant,  the 
messages  enjoy  exponential  service  times,  and  if  we  wish  to  generate  other, 
say  Erlang,  service  times  the  device  of  phases,  or  extra  artificial  compart- 
ments, may  be  used,  as  in  Gaver  and  Lehoczky  (1976). 

If  a message  in  progress  on  a channel  is  interrupted  by  the  arrival  of 
another  message,  both  are  assumed  instantly  destroyed,  and  the  message  initi- 
ators are  transferred  to  a re- try  population,  R . A message  that  has  entered 
R changes  its  status  at  time  t with  probability  v(t)dt  + o(dt)  : with 
probability  a(t)  the  change  of  status  implies  an  attempt  to  re-transmit 
on  a channel,  and  with  probability  a(t)  = l-a(t)  the  change  of  status 
implies  loss  --  the  message  may  no  longer  be  worth  transmitting. 

Our  representation  of  the  above  setup  is  in  terms  of  the  following  state 
variables. 

Q(t)  = the  number  of  messages  being  transmitted,  i.e.  occupying  channels, 

at  time  t ; clearly  0 < ^(t)  ^ c . (2.1) 

JR(t)  = the  number  of  messages  in  R at  time  t ; 0 < £(t)  < 00  . 

L(t)  = the  total  number  of  messages  that  have  been  lost  by  time  t ; 

0 < L(t)  < oo  . 

The  state  of  the  system  is  thus  (£(t),  R(t),  lit))  , a discrete  vector- 
valued Markov  process,  according  to  our  earlier  assumptions,  he  shall 
characterize  this  process  approximately  when  c-*°°  , and  shall  in  particular 
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treat  £ , £ , and  £ as  continuous  stochastic  processes,  in  fact  as 
diffusion  processes,  see  Arnold  (1974). 

Here  is  a formalization  of  the  transitions  described  earlier: 


Transition 

Probability 

(Tn 
R,  L) 

t,  t+dt) 

+ (fitl . £.  L) 

a(t) 

cX(t)[l-  — ]dt 

* <£-’•  £t2'  0 

a(t) 

cX(t)  c dt 

+ (£-'  • £•  i.) 

y(t)£(t)dt 

♦ (a+i • - l) 

Q 

v(t)a(t)[l -^]R(t)dt 
c ^ 

+ (fi-l . R+l . L) 

£ 

v(t)a(t)^£(t)dt 

+ (2,  R-l.  l+l ) 

v(t)a(t)£(t)dt 

(2.2) 


Other  transitions  have  negligible  probability  of  occurring  in  (t,  t+dt)  . 
Now  in  principle  one  can  write  down  Kolmogorov  equations  for  the  transition 
probabilities  of  the  (£,  R,  J.)  process  and  solve  them.  However,  such  a 
solution  must  inevitably  be  numerical.  Here  we  shall  write  down  an  approxi- 
mate system  of  Ito  stochastic  differential  equations,  see  Arnold  (1974).  to 
describe  process  evolution,  and  from  them  deduce  certain  useful  information 
valid  when  c is  large.  We  write,  on  the  basis  of  (2.2), 

Q(t)  Q(t)  Q(t) 

dQ(t)  = (cX(t)[l-'^ — -]  - cA(t)^ y(t)£(t)  + v(t)a(t)[l — ]R(t) 

£(t) 

- v(t)a(t)-£— R(t)}dt 


* 

* 


Qjt) 


+ /cX(t)[l-~— ] dW -| ( t ) - /cA(t)—  dWg(t) 

/ / 5TTJ 

" /n(t)£(t)  d W 3 ( t ) + J v(t)a(t)[l - — — ]R(t)  dW^t) 

/ srn — 

- /v(t)a(t)— R(t)  dWg(t)  , 


(2.3) 


hi 


■ w ■> 


w * '*.*•  . .il 1 


V fv  V 


rr  * 


Q(t)  Q(t) 

dR(t)  = {2cX(t)=- v( t)a( t) [1  - ]R(t) 

^ C C rsj 

Q(t) 

+ v(t)a(t)r::- — R(t)  - v(t)a(t)R(t)}dt 

C m 

/ oltT  / QTt) 

+ 2/cX(t)^—  dW2(t)  - /v(t)a(t)[l-^— ]R(t)  dW^t) 

/ QTt)  / 

+ /v(t)a(t)^- R(t)  dW5(t)  - /v(t)a(t)R(t)  dWg(t)  , 


dL(t)  = v(t)a(t)R(t)dt  + /v(t)a(t)R(t)  dW,(t) 


where  (W^(t),  ost,  i=l,  2,  ..6)  is  a 6 dimensional  standard  Wiener  process 
whose  components  are  independent. 

The  rationale  for  writing  (2.3)  is  as  follows.  Each  term  in  the  brack- 

etted  part  of  the  expression  for,  say,  dQ  , is  a component  of  the  drift 

om 

or  expected  change  in  Q between  t and  t+dt  ; e.g.  cX(t)[l — 1 L]  is 

the  expected  increase  in  Q caused  by  a newly  arrived  message  immediately 

Q ( t ) 

reaching  an  empty  channel  (line  1 of  (2.2)),  while  cX(t)^ — represents 

the  decrease  in  Q caused  by  a new  message  that  arrives  and  chooses  a busy 
channel,  only  to  be  immediately  transferred  to  V , along  with  the  message 
in  progress  (see  line  2 of  (2.2)).  The  coefficients  of  the  Wiener  process 
terms  are  proportional  to  the  standard  deviations  of  the  motions  in 
(t,  t+dt)  corresponding  to  the  drift  terms. 

Next,  introduce  an  expansion  of  (Q,  R,  L)  to  terms  of  order  Acf  . 

m ***  m 

Define  the  stochastic  noise  processes  by 
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Q(t)  - cq(t) 
X(t)  = 

A 


Y(t)  = 


R(t)  - cr(t) 


(2.4) 


L(t)  - cJi(t) 


Z(t)  = 


the  vector  (q(t),  r(t),  &(t))  = lim  (£(t)/c,  JR(t)/c,  JL ( t )/c ) ; the  exis- 

C-KO 

tence  of  these  limiting  functions  is  suggested  by  results  of  Kurtz  (1971 ) . 

An  application  of  Ito's  lemma,  Arnold  (1974),  p.  90,  provides  a sto- 
chastic differential  equation  (s.d.e.)  for  (£,  1)  as  follows  (hereafter 

we  do  not  explicitly  express  the  t-dependence) : 

( dX  ] (-  (2A  + y + 2avr)  av(l-2q)  ] [ X ] 


2A  + 2avr 


-av(l-2q)  - av 


Y + 


q'  - A(l-2q)  + yq  - avr(l-2q) 

/c  r*  - 2Aq  + avr(l-2q)  - avr  + 


(2.5) 


Si ' - avr 


/A(l-q)  -Acf  -/yq^  Avr 


-/avrq  0 


2/\iq  0 /av(l -c 


/vaqr  -/var 


0 0 0 


0 /var 


Now  the  coefficient  of  /c  must  be  identically  zero,  for  otherwise  our 
system  of  equations  does  not  converge.  This  leads  to  the 
Deterministic  Equations: 

q ’ ( t)  = A(l-2q)  - pq  + avr(l-2q) 

r'(t)  = 2Aq  - avr(l-2q)  - avr  (2.6) 

8. 1 (t)  = avr 

The  solutions  to  these  equations,  which  must  be  obtained  by  numerical  inte- 
gration, provide  a deterministic  approximation  to  (Q , R,  L)  --  in  fact, 

^ 

(cq,  cr,  cl)  Z ( E[Q] , E[R],  E[L"|) . 

^ 


Stochastic  Equations: 

In  view  of  (2.6),  (2.5)  now  appears  in  the  form 


dX 

X 

dY 

= \ 

Y 

** 

dZ 

1 

. ~ > 

dt  + B.dW 


(2.7) 


A*  (a  3 x 2 matrix)  being  identified  as  the  coefficient  of  (X,  Y,  Z)1  , 
and  (a  3 x 6 matrix)  identified  as  the  coefficient  of  the  Wiener  process 

term,  in  (2.5).  Next  note  that  the  s.d.e.  has  a special  form  if  the  Z 


term  is  omitted;  it  is  reasonable  to  focus  on  X and  Y , for  clearly 

l( t)-*=°  as  t-*»  , since  L represents  cumulative  losses.  Thus  we  shall  be 

interested  in  the  first  two  equations  in  (2.6),  and  in  (2.7)  which  we  now 
write  as 


dt  + D.dW 


(2.8) 
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and  C*  (2  x 2)  is  A.  without  its  last  row,  while  D.  (2  x 6)  is  B. 

without  i ts  last  row.  Now  the  bivariate  stochastic  process  described  by 
(2.8)  is  nonstationary  Ornstein  Uhlenbeck,  since  £t  has  eigenvalues  with 
non-negative  real  parts;  see  Arnold  (1974),  Sec.  8.3,  for  an  account  of  the 
scalar  case,  or  see  Schach  and  McNeil  (1973).  Much  that  is  useful  is  known 
about  this  process.  In  particular,  if  (£(0),  £(0))  = 0 , and  q(0),  r(0) 
are  given,  then  for  all  t>0  (X_(t),  ^(t))  has  a bivariate  normal  distri- 

bution with  mean  0 and  covariance  matrix  E.  which  satisfies  the  differ- 
ential  equation 

L = c.i*  + e.c;  + d.d:  (2.9) 

~t  ^t^t  ~t  t ~t~t 
see  Arnold  (1974),  Sec.  8.2. 

Combining  these  facts  leads  to  the 

Result:  (£,  £)  is  approximately  bivariate  normal  (Gaussian),  as  c-*»  : 

(&(t),  R(t))  ~N(c(q,r),  c£t)  (2.10) 

From  this  expression  it  is  possible  to  estimate  the  probability  that  at 
least  a specified  number  of  channels  are  occupied  at  time  t , and  also  to 
estimate  the  probability  that  there  are  no  more  than  any  specified  number  of 
customers  awaiting  re-try.  Such  quantities  are  useful  measures  of  system 
performance.  Of  course  both  the  deterministic  equations  (2.6)  and  the 
equation  (2.9)  for  the  covariance  matrix  £t  must  be  solved  numerically, 
but  this  should  be  a less  difficult  step  than  is  the  simulation  of  such  a 
sys  tern . 

Steady-State  Results: 

Suppose  A(t)->-A  , p(t)-*y  , v(t)-n>  , a(t)-Kx  all  positive  constants, 

as  t-«°  . In  this  case  it  may  happen  that  q(t)-*q  and  r(t)-»-r  , the  latter 
satisfying  (2.6)  with  derivatives  set  equal  to  zero: 
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0 = A(l-2q)  - yq  + avr(l-2q) 
0 = 2Aq  - avr(l-2q)  - avr 


(2.11) 

(2.12) 


* 

* 


and  these  are  immediately  solved  to  give 


= (l+2p)  - /(l+2p)2-  8ap  = 


2p 


4a 


(l+2p)  + /(l+2p)2-  8pa 


r = ^ 
a n 


(2.13) 


where  p = A/y  and  n = v/y  , and  provided  that  the  numerical  results  ob- 


tained are  feasible:  0<q<l  , 0<r«»  . From  (2.10),  it  is  necessary  that 


1 


q < 2 in  order  for  a steady  state  to  exist. 


When  steady  state  conditions  hold  Zt=0  , so  the  steady  state  co- 
variance  matrix,  Z , is  the  unique  positive  definite  solution  of 


-DD 1 = CZ  + ZC1 


where 


C = 


- ( 2A  + p + 2avr) 
2A  + 2avr 


av(l-2q) 
-av(l-2q)  - av 


and 


DD*  = 


A + yq  + avr 
-(2Aq  + avr) 


-(2Aq  + avr) 
4Aq  + vr 


writing 


Z = 


a a 
11  12 


a a 
12  22 


we  may  express  the  solution  of  (2.14)  as  follows: 
fb(a+b-d)  2b2  b2 

bd  2ab  ab 


' ' 

a 

1 1 

- A"1 

a 

1 2 

a 

22, 

2ad  a(a+b)  - bd 


A + yq  + avr 
-(2Aq  + avr) 


4Aq  + vr 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


j 


where 


a = 2A  + y + 2avr  , d = 2A  + 2avr 

(2.19 

b = av(l-2q)  , A = 2ab(d-a-b)  + 2b2d  . 

The  values  of  q and  r are  available  from  ( 2 . 1 J y . From  these  results 
actual  numerical  expressions  for  the  steady  state  probability  distributions 
of  Q and  R may  be  computed.  Notice  that  the  long-run  loss  or  defection 
rate  of  the  final  equation  of  (2.6)  i.e.  cil 1 (°° ) = cavr  , where  r comes 
from  (2.13). 


3.  Model  2:  Intelligent  Re-tries 

Many  variations  are  possible  on  the  process  that  leads  to  Model  1. 

These  are  of  interest  because  they  represent  possibilities  for  improved 
performance,  or  because  they  describe  real  system  behavior. 

In  the  model  of  the  present  section  we  suppose  that  re-tries  are  in- 
telligent: a message  in  category  R that  attempts  to  seize  a channel 

again  does  not  do  so  at  random,  with  a chance  of  destroying  a message  in 
progress.  To  represent  this  change,  remove  in  (2.2)  the  possibility  of  the 
transitions  involving  a(t):  (Q,  R,  L)  ■>  (Q+l , R-l , L),  -*■  (Q-l , R+l , L), 

and  -*•  (Q,  R-l,  L+l ) , replacing  with  the  transition  rate 

Q 

Q(t),  R ->  Q(t)  + 1,  R(t)  - 1 , probability  v(t)R(t)[l  - -]dt  (3.1) 

Effectively  this  change  in  Model  1 prevents  losses,  and,  when  a re-try 
prepares  to  access  an  occupied  channel,  it  senses  the  presence  of  the  on- 
going message,  and  refrains  before  destruction. 

This  model  can  be  analyzed  by  the  technique  described  for  Model  1, 
namely  that  of  setting  up  approximating  Ito-type  stochastic  differential 
equations,  and  then  carrying  out  an  expansion  for  c -*■  00  . The  results  are 
as  follows. 

Deterministic  Component: 

We  find  that  for  Model  2 

q'(t)  = A(l-2q)  - pq  + vr(l-q) 


r'  (t)  = 2Aq  - vr(l-q) 


Stochastic  Component: 


Use  of  the  representation  (2.4)  together  with  the  stochastic  differ- 
ential equations  and  Ito's  lemma  yields 


+ hd"t  ’ 


- *;\.v.vrv  * 


the  s.d.e.  for  a non-stationary  Ornstein-Uhlenbeck  process.  In  this  case 
(suppressing  t-dependence) , 


5t" 


-(2A  + y + 

vr) 

v(l-q) 

2A  + vr 

-v(l-q) 

*^TT-qT 

-Aq 

-A? 

0 

2Aq 

0 

/vr(l-q) 

-/vr(I-q) 


(3.4) 


(3.5) 


and 


Wt  = (W-,  (t) , W2(t),  W3(t),  W4(t))'  , (3.6) 

a 4-dimensional  standard  Wiener  process  with  independent  components. 

Once  again  the  equation  (2.9)  may  be  solved  for  the  covariance  function 

, with  At  substituted  for  Ct>  and  in  place  of  Dt  . And  once 

again  we  state  the 

Result:  (Q,  R)  is  approximately  bivariate  normal  as  c -*■  °°  : 

(Q,  R)  ~ W(c(q,r),  clt)  (3.7) 

Steady-State  Results: 

If  A(t)  -*•  \ , y(t)  -*•  y , v(t)  -*  v as  t -*•  <*>  , then  it  may  happen  that 

q(t)  •*  q , and  r(t.)  ■+  r , these  satisfying  the  equations  (3.2)  with 

q1  = r'  = 0 . The  solutions  are  easily  obtained,  and  are 


n_A_  _ _ 2A2  . 2p2  . 

q y ~ P ’ v(y-X)  _ n(l-p)  * 


(3.8) 


apparently  it  is  necessary  that  q = p < 1 for  steady  state  to  occur. 

Note  that  in  Model  1 it  was  necessary  that  q < 1/2  in  order  that  steady 
state  conditions  occur.  If  we  compare  Models  1 and  2 when  a = 1 we  find 
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Model  1 Model  2 

q = p » (o  < p < ^)  q = p»  (o  < p < l ) 

_ _ 2p2  y.  = 2p2 

n(1-2p)  n(l-p) 

Clearly  the  setup  of  Model  2 results  in  considerably  smaller  expected  wait 
per  unit  time  (approximated  by  c(q+r))  than  does  that  of  Model  1.  The 
improved  service  must  be  purchased  in  return  for  investment  in  the  busy 
channel  sensing  capacity.  The  Model  2 can  be  compared  to  the  model  of  Gaver 
and  Lehoczky  (1976). 

The  steady  state  covariance  matrix,  Z , is  obtained  from  (2.14), 
replacing  C by  A , and  D by  B . We  find  that 


a 

i i 

b(a+b-d) 

2b2 

b2 

A + yq  + vr(l-q) 

0 

- A'1 

bd 

2ab 

ab 

-(2Aq  + vr(l-q)  ) 

1 2 

a 

d2 

2ad 

a(a+b)-bd 

4Aq  + vr(l-q) 

22 

. 

. - 

where 


a = 2A  + y + vr  , b = v(l-q) 
d = 2A  + vr  , A = 2ab(d-a-b)  + 2b2d 


4.  Model  3:  Transitory  Version  of  Model  1. 


Consider  next  a transitory  service  system  version  of  the  basic  Model  1. 
There  are  N messages  to  be  transmitted,  and  each  initially  is  transmitted 
independently  and  at  a time  having  absolutely  continuous  distribution  function 
F(t)  , with  F(0+)  = 0 , and  density  f(t)  . Once  a given  message  is  trans- 
mitted, no  more  appear  from  its  particular  source,  so  even  though  some 
messages  enter  the  re-try  population  one  or  more  times,  the  traffic  gradu- 
ally fades  away;  the  problem  is  fundamentally  non-stationary . For  analysis 
of  a similar  problem  see  Gaver,  Lehoczky,  and  Perlas  (1975). 

The  following  state  variables  are  required. 

I ( t)  = the  number  of  message  arrivals  that  have  occurred  by  time  t ; 

0 ^ I(t)  < N . 

Q(t)  = the  number  of  messages  being  transmitted  at  time  t ; 

0 < Q(t)  s c . 

R(t)  = the  number  of  messages  in  the  re-try  population  at  time  t . 

We  assume, as  we  did  in  formulating  Model  1,  that  if  a newly  arriving  message 
encounters  a channel  --  selected  at  random  — held  by  a message  in  progress, 
then  both  are  instantly  "destroyed,"  and  immediately  join  R . No  messages 
are  ever  lost.  Messages  in  R re-try  at  rate  v(t)R(t)  . 

The  transition  scheme  is  given  below.  The  individual  message  arrival 
rate  at  time  t is  seen  to  be  X(t)  = f ( t )/ [1 -F( t ) ] . 


Transition 
(t  to  t+dt) 

(I,  Q,  R)  - (1+1,  Q+l , R) 

- (1+1,  Q-l,  R+2) 

- (I,  Q-l,  R) 

- (I,  Q+l,  R-l ) 

- (I,  Q-l,  R+l ) 


Probability 


A(t)[N  - I(t)][l  - £-]dt 

Q 

A(t)[N  - I(t)]^dt 
y(t)Qdt 

Q 

v(t)R[l  - T-]dt 

~Q  C 

v(t)R^dt  . 


(4.1) 
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r'r*wr7=rr- 


. >■  > . 


» » 


T 


From  the  above  we  can  write  down  the  approximating  stochastic  differential 
equations  analogous  to  (2.3  ).  These  are,  suppressing  t-dependence  state- 
ments. 


dl  = X [N  - I]dt  + A[N  - I][l  - dW-j  + A[N  - I]^  dW2  (4.2) 


Q Q Q Q 

dQ  = U[N  - I][l  - £■]  - A[N  - !]§■  - yQ  - vjg-  + vR[l  - ^]}dt 


/ Q / Q _ 

+ A[N  - I][l  - §•]  dW]  - A[N  - I]^r  dW2  " d~3 


fQ  / 9 

- /vR-  c'u^  + /vR[l  - -]  dW5  , 


Q Q Q 

dR  = {2A[N  - IE-  + vR=-  - vR[l-  ^]}dt  + 

~ ~ L ■VL  ~ ^ 


Once  again  we  study  the  behavior  of  the  system  as  system  parameters,  in  this 
case  both  N , the  initial  number  of  messages,  and  c , the  number  of 
channels  become  large.  In  fact,  we  relate  these  parameters  as  follows: 
c = BN,  B being  a positive  constant.  Introduce  the  noise  processes 


X(t)  = 


I(t)  - Ni(t) 


Y(t)  = 


Q(t)  - Nq(t) 


(4.3) 


Z(t)  = 


R(t)  - Nr(t) 


where  (i(t),  q(t),  r(t))  = lim  (I(t)/c,  Q(t)/c,  R(t)/c)  . 

N-x» 

Then  an  application  of  Ito's  lemma  and  identification  of  terms  of  order 
i/N  and  N yields  the  deterministic  and  stochastic  components  of  the  process, 


Deterministic  Equations: 


i 1 (t)  = X(1  - i) 

q ' ( t)  = A(1  - i)(l  - q/B)  - A(1  - i )q/B  - yq  - vrq/B  + vr(l  - q/B)  . 


r'(t)  = 2A(1  - i)q/B  + vrq/B  - vr(l  - q/B) 


(4.4) 


From  the  first  equation  and  the  definition  of  X it  follows  immediately  that 
A(l-i)  = f , the  density  of  the  arrival  distribution.  Note  that  when  this 
substitution  is  made  in  the  second  and  third  equations  the  latter  are  precisely 
of  the  form  of  the  corresponding  equations  of  Model  1: 


Correspondence  Between  Model  1 and  Model  3 
Arrivals  and  Re-tries 


Model  1 


Model  3 


f(t)/B 

v(t)/B 


In  the  present  model  f(°°)  = 0 , and  hence  the  arrival  rate  is  eventually 
zero,  and  some  time  afterwards  the  system  is  completely  empty.  The  present 
model  possesses  no  steady  state  solution,  and  to  learn  about  system  status  at 
various  times  it  is  necessary  to  solve  (4.4)  numerically. 


Stochastic  Equations: 

The  expansion  technique  shows  that  the  noise  process  satisfies 


where 


dX 

■i 

X 

dY 

= At 

Y 

dZ 

Z 

+ 5 A 


-AO  0 

-A(l-2q/B)  -2f/B  + y+  2vr/B  v(l-2q/B) 

-2Aq/B  2f/B  + 2vr/8  -v(l-2q/B) 


and 


✓fU-q/B) 

/fq73 

0 

0 

0 

II 

4-> 
CO  i 

/f(1-q/B) 

-/Fq7B 

-4Iq 

-/vrq/B 

/vr(l-q/8) 

0 

2/fq7B 

0 

/vrq/8 

-/vr(l-q/8) 

(4.4) 


(4.5) 


(4.6) 


and  Wt  is  a 5-dimensional  standard  Wiener  process  with  independent 
components. 


Result;  (I,  Q,  R)  is  approximately  normal  as  N (hence  c)  -+•«*>  : 

(I,  Q.  R)  ~ W(  c(f,  q,  r),  cEt  ) 

Here  f t satisfies  the  differential  equation 


E*  = A.E.  + E.A!  + B.B; 
~t  ~t~t  ~t  t ~t  t 


and 


5t!t 


(4.7) 


2fq/B 


f f(l-2q/B) 

f(l-2q/B)  f(l-2q/B)  + yq  + vr  -(2fq/S  + vr) 

2fq/B  - ( 2fq/3  + vr)  4fq/B  + vr 


(4.8) 
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As  might  be  anticipated,  the  fact  that 

Var[I(t)]  =o  (t)  + F(t)[l-F(t)] 
~ 1 1 

may  be  deduced  from  (4.7). 
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Suppose  now,  as  may  be  quite  realistic,  that  a message  attempting  to 
transmit  on  a channel  does  so  until  completion  before  it  discovers  that  it 
has  been  "destroyed,"  i.e.  garbled  by  others--which  it  also  contributes  to 
garbling.  The  analysis  at  once  becomes  much  more  complex  because  a channel 
may  contain  any  number  of  destroyed  messages.  Our  formulation  is  that  of 
Model  1,  save  for  the  change  described  above. 

The  following  state  variables  describe  the  system. 

Q(t)  = the  number  of  channels  carrying  good,  i.e.  ungarbled  or  unde- 
stroyed messages  at  time  t . 

Sk(t)  = the  number  of  channels  carrying  exactly  k messages  that  are 
destroyed  at  time  t ; clearly  SQ(t)  e Q(t)  , and 

0 < Q + ? S.  < C. 

~ k=rK 

R(t)  = the  number  of  messages  in  R at  time  t . 


The  transition  probability  scheme  is  summarized  next.  We  write 


(g.  §1 


Transitions 


S2  , • " ' 

..  s 

k ?> 

-► 

(Q+l, 

h* 

• sk , . 

...  R) 

(Q+i. 

§1 » 

...»  St  , . 

...  R-l ) 

(Q-l. 

....  sk , . 

...  R) 

(Q-l. 

~r 

$2+1 . 

. sk,.... 

R-l) 

(Q-l. 

§1 » 

$2+1 , . . . 

* 

R) 

(9*  §1 

i»** 

. , sk+i , . 

R-l) 

for 

k = 

1,  2,  3, 

• • • 

(9* 

i * * ■ 

sk-i. 

5k+l  • — • 

R) 

(9*  9- 

r1- 

....  sk>. 

R+l ) 

(Q » S-| , . . . , Sk+1 » 5k+-j  » 


Probabilities 
Q+S+ 

c A [ 1 - — — ]dt 
Q+s+ 

vR[l  — - — ]dt 

uqdt 
Q 

vR-^dt 
~c 

Q 

cAF 

vR^dt 
~c 

Sk 

cXrdt 

~i 

y — dt 
c 

R+l)  y(k+l)Sk+1 


(5.1) 


There  will  be  a denumerable  infinity  of  such  transitions.  While  it  is  in 
principle  possible  to  carry  out  the  expansion  technique,  we  shall  content  our- 
selves with  a brief  discussion  of  the  deterministic  equations  analogous  to 
(2.6).  Following  the  example  of  Section  2,  we  can  write  down  these  differ- 
ential equations  for  the  limits  as  c-*»  of  Q/c  , R/c  , and  S^/c  (k=l,2,...), 
denoted  by  q , r , and  sk  : 
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V.*  .V 

■ - - - - — - - 


q'  = - (X  + vr)(l  - s+  - q)  - (y  + X + vr)q 

00 

r'  = - vr  + y £ ks. 

k=l  K (5.2) 

sk  = y(k+l)sk+]  + (X+vr)sk_i  - (k+X+vr)sk  , k=l ,3,4,5, .. . 

S2  = 3ySg  + (X+vr) (s-j+q)  - (2y+X+vr)$2 

The  steady-state  values,  which  exist  under  circumstances  to  be  derived,  are 
obtained  by  solving  the  above  system  of  equations  with  the  derivatives  set 
equal  to  zero.  To  simplify,  first  divide  through  by  y , putting 
p = X/y  , n = v/y  . Af  3r  sortie  algebra  it  is  found  that 

q = P (5.3) 

p + nr  = 0 , so  r = (e  - p)n"^  (5.4) 

where  0 is  the  smallest  solution  of  the  equation 


x -x 
e = p 


TTx 


(5.5) 


provided  one  exists. 
s1  = p0 


sk  = 


e-00k 
k! 


(5.6) 


k > 2 


Further  analysis  shows  that  in  order  for  a steady-state  to  exist, 

p < -- — ■ -■  e'^  ' ])/2  = 0.20588... 

+ 1 

Notice  that  this  value  is  much  smaller  than  the  just-intolerable  value  of 
1/2  derived  from  Model  1.  Clearly,  transmission  to  completion  of  messages 
provides  opportunity  for  many  more  transmissions  to  be  destroyed. 
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6.  Model  5 • Transmission  to  Completion  with  Intelligent  Re-tries 

A natural  variation  of  the  previous  model  is  obtained  by  insisting 


on  transmission  to  completion,  but  also  allowing  for  an  intelligent  re-try 
capability,  as  in  Model  2.  This  means  that  the  transitions 
-«H.  S i » $2+1 , . • t SR , . . , R-l ) , and  "^(Q,  S-j,...,  R-l)  are 

ruled  out,  i.e.  have  probability  zero  in  this  model.  Consequently  the  tran- 
sition rates  of  (5.1)  apply,  with  this  change.  Confining  attention  to  the 
deterministic  differential  equations  it  may  be  shown  that  these  have  the 
form  given  below. 

q'  = (A+vr) (1  - S+  - q)  - (A+y)q 


r*  = -vr(l  - s - q)  + y £ ks. 

k=l  K 


These  are 


q = P 

r = 1 gLU.+p)eP-'l]. 

n l-(l+p)(ep-l ) 


sr  p 

k 

sk=  (1+P)kT  * k=2,3, . . . 


(6.1) 


sR  = (k+1)  ysR+-|  + As k_ i - (A+ky)sR 
s^  - Sys^  + As-j  - (A+2y)$2  + Aq  . 

Once  again  let  p=A/y,  n=v/y,  and  solve  for  the  steady-state  values. 


(6.2) 

(6.3) 

(6.4) 

(6.5) 


The  condition  for  existence  of  a steady  state  is  that  the  denominator  of 
(6.3)  be  positive,  which  translates  into  the  requirement  that  p < 0.50855...  . 
The  latter  value  may  be  contrasted  to  the  just-tolerable  value  of  unity 
derived  for  Model  2. 
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7.  Numerical  Results  for  Model  1. 


In  order  to  check  the  quality  of  the  diffusion  approximations,  a 
simulation  program  was  written  that  realizes  the  two-dimensional  Markov 
process  describing  Model  1.  The  latter  was  then  exercised  through 
5 x 10^  state  changes  for  several  values  of  the  offered  load.  A/p  , 
and  for  several  channel  numbers;  the  parameter  a = 1 , and  v/p  = 0.08 
in  this  experiment.  The  results  were  generally  supportive  of  the  approxima 
tion,  as  is  seen  by  examining  the  following  tables.  Note,  however,  that 
assessments  of  the  mean  and  variance  of  the  number  of  active  channels,  Q , 
are  generally  in  closer  agreement  than  are  those  for  the  size,  R , of  the 
re-try  population. 


TABLE  1. 


Means  and  Variance-Covariance  Values  by 
Simulation  and  Diffusion  Approximation 


c = 1 0 ; v/y  =2.0 


A/y:  0.48 

0.40 

0.32 

0.20 

E[Q]  simulat. 

4.79 

3.99 

3.19 

1.99 

diffus. 

4.80 

4.0 

3.2 

2.0 

Var[Q] 

2.52 

2.52 

2.38 

1 .76 

2.50 

2.52 

2.39 

1.77 

E[R] 

65.13 

9.31 

3.29 

0.75 

57.6 

8.0 

2.84 

0.67 

Var[R] 

1470. 

50.26 

10.62 

1.57 

1466. 

42.81 

8.80 

1.35 

Cov[R,Q] 

0.034 

0.124 

0.166 

0.149 

0.038 

0.143 

0.177 

0.150 

« 

* 


TABLE  2. 


Means  and  Variance-Covariance  Values  by 
Simulation  and  Diffusion  Approximation 

c = 25  ; v/y  = 2.0 
X/y:  0.48  0.40 


E[Q]  simulat. 


11.95 


9.97 


8.  Summary  and  Conclusions 


In  this  paper  it  has  been  shown  that  an  approximate  approach,  using 
stochastic  differential  equations,  is  effective  for  modeling  an  element  of 
a complex  communication  system.  The  approximation  improves  as  c , the 
number  of  channels  available  for  transmission,  becomes  large.  The  adequacy 
of  the  model  under  such  conditions  is  suggested  by  the  theoretical  results 
of  Barbour  (1974)  and  Kurtz  (1971);  numerical  solutions  of  selected  systems 
also  attest  to  the  adequacy  of  the  approximation.  The  authors  wish  to 
gratefully  acknowledge  the  research  support  of  the  Office  of  Naval  Research. 
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